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Dipole excitations in neutron-rich nuclei below the neutron threshold are investigated. 
^ _ The method is based on Hartree-Fock-Bogoliubov (HFB) and Quasiparticle-Phonon Model 
O ' (QPM) theory. Of our special interest are the properties of the low-lying 1~ Pygmy Res- 
^ ■ onance and the two-phonon quadrupole-octupole 1~ states in Sn-isotopes including ex- 
■ ploratory investigations for the experimentally unknown mass regions. In particular we 
^ ' investigate the evolution of the dipole strength function with the neutron excess. The use 
of HFB mean-field potentials and s.p. energies is found to provide a reliable extrapolation 
into the region off stability. 

Ch ■ 1- Introduction 
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The understanding of the structure of atomic nuclei at extreme isospin is a challenge for 
nuclear physics which has become achievable on a large scale after the modern experimen- 
ts \ tal nuclear structure facilities [1^ came into operation. For the theoretical studies of such 
systems of importance is the search for suitable observables. One promising choice for 
the investigation for example of neutron-rich nuclei is to look after their excitation modes 
and trying to find correlations e.g. of transition probabilities and excitation energies with 
the size and shape of the neutron skin. The appearance of low-energy electric dipole 
strength is a genuine feature of neutron-rich nuclei, seen recently in high-precision photon 
scattering experiments already in stable nuclei with small |2|3j and moderate |3] neutron 
excess. These so-called Pygmy Dipole Resonances (PDR) are located close to the neutron 
threshold, forming in presently accessible medium- and heavy-mass nuclei a clustering of 
states in the region of E^, ~ 5.5-=-8 MeV. Although carrying only a small fraction of the 
full dipole strength these states are of particular interest because they are refiecting the 
motion of the neutron skin against the core of normal nuclear matter. As discussed in 
|3] , their nature is completely different from the representation where proton and neutron 
fluids as a whole move against each other, giving rise to the prominent isovector Giant 
Dipole Resonance (GDR). 

The PDR mode has to be distinguished from the other known low-energy isoscalar 
dipole excitation, namely the two-phonon 1^ states resulting from the anharmonic inter- 
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actions of the lowest 2"*" and 3~ states in a nucleus. A comprehensive collection of data 
on the latter and overview of the widely used experimental methods is given in [S]. The 
anharmonicities are reflecting the intrinsic fermionic structure of the nuclear phonons 
thus deviating from ideal bosons. A model taking into account that nuclear phonons 
are formed by exciting at least two quasiparticles (corresponding to a one particle one 
hole state (Iplh)) which interact by a residual two-body NN-interaction is the QPM |6|7j 
worked out and applied systematically |8|9|lUj . e.g. also in the recent PDR investigations 
in the stable ^^^'^^''Sn|ilj and ^osp]^ gj nuclei. Other applications of the QPM to low- 
energy dipole strength |12|13|14j have led to very good description of data thus giving 
confidence on the reliability of the model for such investigations. 

Most likely, the QPM approach is at present the only available method for a unified 
description of low-energy single and multiple phonon states. The use of schematic interac- 
tions of separable form and with empirical coupling constants is a detail which in principle 
could be improved on. Recently, finite rank approximation for Skyrme interaction was 
incorporated in the formalism of QPM Work in this direction is in progress but 

since the investigations in this paper are mainly directed to an exploratory study of the 
evolution of the PDR modes with neutron excess we use the standard form of the QPM. 

Here, the mean-field part, however, is treated microscopically by incorporating HFB 
results on single particle energies and potentials as input for the QPM calculations. 

The interest on the region around the N=82 shell closure has arisen recently |16|17|18] 
and will further increase with the new experimental opportunities on REX-ISOLDE, 
MAFF and also the coming-up GSI project. In this manuscript we have presented calcu- 
lations for the neutron-rich Sn isotopes with neutron number N=70-=-82 as an exploratory 
investigation in experimentally less or even unknown regions. 

2. Description of the method 

For the extrapolation of the QRPA and QPM calculations into unknown mass regions 
a rehable description of the ground state properties is necessary. Because of the nu- 
merical constraints set by the QPM a semi-microscopic approach is chosen. The model 
Hamiltonian is written as: 



where Hmf = Hgp + Hpair represents the mean-field part describing the motion of inde- 
pendent quasiparticles in a static potential and interacting in the particle-particle (p-p) 
channel through the monopole pairing interaction Hpair- Hmf is to be identified with the 
HFB Hamiltonian discussed in |19|2Uj . For practical purposes, however, we use here a 
re-parameterization of the HFB mean- field in terms of Wood-Saxon potential- Uy/s fitted 
to the HFB mean-field such that the single-particle energies are reproduced. Also, the 
pairing part is simplified by using the conventional constant matrix element approach 
which close to the Fermi level is a reliable approximation. 

H^M -^5M residual interactions taken as a sum of isoscalar and isovector sep- 
arable multipole and spin-multipole interactions in the particle-hole channel. H^ is the 
sum of the multipole pairing interactions in the particle-particle (p-p) channel. 

In practice, for a given nucleus of mass A the depth of the central and spin-orbit 
potentials, radius and diffusivity parameters of Uws are adjusted separately for protons 



H = Hmf + HI'} + HI\, + H^ 
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and neutrons to the corresponding single particle separation energies, the total binding 
energy [7T], the charge radii and (relative) differences of proton and neutron root-mean- 
square (RMS) radii, 

Sr = ^J<r^ >„ - ^J<r^ >p (2) 

indicating the thickness of the neutron skin. The results from this procedure are compared 
to experimental values from the compilation of Audi and Wapstra j2I] and the recent 
investigations of neutron skins by Krasznahorkay et al. by charge exchange reactions [22] , 
respectively. For both observables the agreement of theory and data is very satisfactory 
as it is seen from Fig.l. 
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Figure 1. Ground state properties of the Sn isotopes. In the upper panel the total 
nuclear binding energies per particle are compared to the values from the Audi- Wapstra 
compilation [21]. In the lower panel the differences of proton and neutron rms radii are 
compared to the experimental values obtained from charge exchange reactions in ref . [22] • 



The model Hamiltonian is diagonalized assuming a spherical 0+ ground state which 
leads to an orthonormal set of wave functions with good total angular momentum JM. 
For even-even nuclei these wave functions are a mixture of one-, two- and three-phonon 
components in the following way [IB] : 

^AJM) = lY.R^WQJM^ + E Px^'^W 
\ i Alii 
A2«2 
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where R,P and T are unknown amplitudes, and u labels the number of the excited states. 

The model basis includes the natural parity 1~, 2+, 3~, 4"*", 5~ phonons only. For 
the description of the structure of each of the excited states we use a wave function 
from Eq. (jS)). For the case of 1~ states one-phonon configurations up to Ex=20 MeV 
are included. This allows us to take into account microscopically the influence of core 
polarization on the low-lying 1~ states without any phenomenologically induced effective 
charges. The two- and three-phonon configurations are truncated up to 4.5 MeV for 
the calculation of the quadrupole-octupole 1~ state for the proper comparison with the 
available NRF data ,12j. For the QPM calculations between 4.5-i-8 MeV the two- and 
three-phonon basis is limited to states up to Ex=8.5 MeV and 8 MeV, respectively. 

3. Results 

The new method we have applied for the determination of the ground state properties 
allows to investigate the evolution of the dipole strength function with the neutron excess. 

Of our particular interest are the one-phonon 1~ states below the neutron threshold with 
excitation energies up to 8 MeV in i20^i26gj^ ^^^^ ^ MeV in i28^i32gj^ respectively. 

In fact, the first three 1~ QRPA states enter in this energy region. They are rather 
well separated from the other higher-lying one-phonon 1~ states by an energy gap of 
more than 1.3 MeV. The first QRPA state contains mainly a neutron two-quasiparticle 
(~ 99.7%) component while the mixing between different two-quasiparticle neutron and 
proton configurations (the contribution of the protons is about 0.2% only) becomes more 
important for the second and the third QRPA states. These states we identify with 
the PDR which in the case of our QRPA calculations is located in the energy region 
5.9^7 MeV. The obtained neutron structure of the 1~ QRPA states is in agreement with 
the explanation of the PDR in other nuclei obtained before by DFT |23|, RQRPA HH] 
and QPMg theory. 

This remarkable stability of the wave functions is similar to what is found in theoretical 
studies of the GDR. In medium and heavy nuclei the model independent Thomas-Reiche- 
Kuhn limit of the dipole EWSR is found to be almost exhausted already by RPA or QRPA 
calculations, indicating the dominance of Iplh or two-quasiparticle structures and agree- 
ing with experiment. However, we emphasize that the PDR states are of predominantly 
isoscalar (or eventually mixed isospin) character as already pointed out in Hence, 
the PDR states cannot be considered as belonging simply to the low-energy tail of the 
isovector GDR. Rather, this dipole modes are of a genuine character which cannot be 
deduced by extrapolations from the GDR region. 

As a general result we find a correlation of the total PDR strength and the number 
of the neutrons (presented on Fig.2a). The value of the total PDR strength increases in 
going to the heavier tin isotopes, up to ^^^Sn. At the same time the centroid energies 
show an opposite behavior (Fig. 2b). 

We have extended the QRPA calculations in terms of QPM, which allows us to inves- 
tigate the properties of the lowest-lying two-phonon quadrupole-octupole 1~ states 
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Figure 2. QRPA calculations on the PDR a)total strength and b)centroid energy as a 
function of the neutron access in ^^''"^^^^Sn isotopes. 



in i20-^i30gj-^ isotopes as well. 

A proper description of the quadrupole J'^=2f and octupole J'^=3{ states is important 
with respect to the two-phonon J'^=l~ state. The calculated energies and transition 
probabilities of the first 2+, 3~ states are presented in Table 1, where also the known 
experimental data are shown A good agreement between the calculations and the 
experiment is obtained. 

The properties of the lowest lying 1~ states have been investigated in the frame of 
QPM. The energy varies from E^.=3.203 MeV to E^.=4.115 MeV in ^^o-iaog^ (g^g r^^^jg 
1). In these isotopes the structure of the lj~ states is predominantly by about 88-^93% 
of two-phonon quadrupole-octupole character. Since with increasing neutron number the 
energy of the latter becomes larger, certain other, higher-lying, two-phonon configurations 
start to contribute as well. The three-phonon configurations are most important for the 
investigated lower mass number Sn isotopes. This effect is connected with the decrease 
of the collectivity of the nuclear excitations when we come to the closed N=82 shell. The 
El transitions from the two-phonon 1~ state to the ground state and between 3^ and 2i 
excited states are an example of 'boson forbidden' transitions with one- or two-phonon 
exchange [HI. The QPM calculations on the B(E1; 0+ 1^) and B(E1; 3^ ^ 2^) 
transition probabilities are in a good agreement with the available experimental data. 

A comparison between the QRPA and QPM dipole strength distribution below the 
neutron threshold in ^^^Sn and ^^°Sn is presented in Fig. 3. The lowest QPM 1~ state 
has no QRPA counterpart because, as discussed above, it is of quadrupole-octupole two- 
phonon structure. While in QRPA the dipole strength in the investigated Sn isotopes is 
concentrated in three dipole states the QPM phonon interactions lead to fragmentation - 
in this case about 80 states - and dissipation seen as a reduction of the total one-phonon 
content of the multi-phonon eigenstates. 

These quenching and fragmentation effects are directly observable in the B(E1) transition 
probabilities to the ground state which are determined by the one-phonon components 
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Table 1 

QPM results for the energies and the reduced B(E1), B(E2) and B(E3) transition prob- 
abilities of the first 1~, 2"*" and 3~ states in ^^°"^^^°Sn isotopes. A comparison with the 
experimental data is presented TP. 



Nucl. Energy Trans. B(E1; J^,') [IQ-^ e^fm^] 
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Figure 3. QRPA and QPM calculations of the El strength distribution in a/^^^Sn and 
b/^^°Sn below the neutron threshold. 



of the excited states: reductions in the one-phonon state amplitudes are reflected in a 
corresponding suppression of the ground state transition matrix elements. In ^^^Sn and 
^^°Sn about 85% of the FDR one-phonon strength is exhausted in the energy intervals 
up to 8 MeV and 7.5 MeV, respectively. Compared to the calculations in ^1] we find 
less fragmentation because the size of our two- and three-phonon configuration spaces are 
somewhat smaller. However, the total El strength in the investigated energy interval is 
in very reasonable agreement with the available experimental data in ^^^Sn. 



4. Summary and Conclusions 

A new method based on HFB and QFM theory was applied for the study of the neutron- 
rich tin isotopes. From our calculations on ^^°Sn ^'^^Sn we obtained low-energy dipole 
strength in the energy region 5.94-6.9 MeV concentrated in a narrow energy interval such 
that a pygmy dipole resonance (FDR) can be identified. The correlation of the FDR 
excitation energy and transition strength with the neutron excess was investigated. Its 
strength increases with the neutron number while the centroid energy decreases. By 
comparison to the radii of HFB ground state densities a correlation with the size of the 
neutron skin could be identified. An important step in understanding the dipole spectra 
is to disentangle the FDR states from the low-energy two-phonon dipole states. In our 
calculations this was achieved by using the QFM approach with up to three-phonon 
configurations. 

For the experimentally unknown nuclei the excitation energy of the two-phonon 1~ 
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state and B(E1) transition probability are predicted. The method will be applied for 
other regions of nuclei. 

The authors of the paper wish to acknowledge the discussions they had with U. Kneissl 
and H.H. Pitz. This work is supported by DFG, contract Le439/5. 
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